Monte Carlo simulations reveal the straightening up of an end-grafted flexible chain 

with a rigid side chain 
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We have studied the conformational properties of a flexible end-grafted chain (length A*') with 
a rigid side chain (length S) by means of Monte Carlo simulations. Depending on the lengths A'^ 
and S and the branching site, b, we observe a considerable straightening of the flexible backbone as 
quantified via the gyration tensor. For b = N, i.e. when attaching the side chain to the free end of 
the flexible backbone, the effect was strongest. 



I. INTRODUCTION 

Polymers anchored on a locally flat substrate are of 
great importance for functionalized surfaces in the ma- 
terial sciences [J and in biology 0- A particular ex- 
ample for the latter is the protective extracellular ma- 
trix of living cells that is composed of the flexible poly- 
mer hyaluronic acid (HA) to which semi-flexible aggrecan 
chains are attached via linker proteins [1, 0] ■ The HA- 
aggrccan has recently received an increased attention as 
it is not only used to equip cells with a protective layer 
but also plays an important role as gliding surface in the 
articular cartilage of synovial joints Understanding 
the material properties of the HA-aggrecan system on the 
nano- and mesoscale are thus of key importance when 
aiming at constructing biomimetic surfaces that mimic 
the function of natm'al cartilages. In fact, on the singlc- 
moleculc scale the HA-aggrecan system can be simplified 
to a flexible polymer, attached to a planar substrate, with 
rigid side chains. 

Flexible, self-avoiding end-grafted polymers without 
side chains have been studied extensively by theory, ex- 
periment and computational approaches [5|,l6|. Typically, 
a large set of end-grafted polymers has been considered, 
where the grafting density cr, i.e. the distance between 
individual chains, was varied ('polymer brush'). In the 
limit of very low cr, the polymer chains can be considered 
as isolated entities, each acquiring approximately a half- 
spherical shape ('mushroom') with a radius comparable 
to the Flory-radius Rp ^ a ■ N^/^ of a coil in a good 
solvent 0. Here, a denotes the size of a monomer. In 
fact, a more detailed description of the polymer shape 
requires the gyration tensor : 
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where r*-'^ is the position of the ith monomer and in- 
dices m, n denote its individual components. Using an 
orthogonal transformation, S can be converted to diago- 
nal form with entries < < denoting the squared 
lengths of the principal axes of gyration while the trace of 
S yields the squared radius of gyration. The associated 



eigenvectors of the gyration ellipsoid are denoted by Vi , 
V2, and V3, respectively. 

The ratios of the principle axes of inertia quantify the 
deviation of the polymer shape from a sphere. Using 
Monte Carlo simulations it was shown that simple ran- 
dom walks (RW) reveal a pronounced asphericity ('as- 
phcricity of the RW) mirrored in the asymptotic ra- 
tios (Ll) : (Ll) : {Lfj ^ 1 : 2.7 : 12.0 [l^. Self- 
avoiding walks (SAW) show an even more aspheric shape 
(1 : 2.98 : 14.0) [lO| that becomes more pronounced when 
attaching the SAW with one end to a flat, solid substrate 
(1 : 3.0 : 14.9 for TV ^ 00) [IH. 

The conformational properties of polymer systems 
comprising side chains branching out from a flexible or 
semi- flexible backbone at varying density ('comb poly- 
mers') have also been considered extensively in theory 
and simulation studies [l2|. The flexibility of the en- 
tire complex has been shown to be strongly influenced 
by the density of flexible side chains with three scal- 
ing regimes in case for a free (i.e. not end-grafted) back- 
bone [1^, [ij] . One of the main results was that attaching 
a large number of flexible side chains stiffens an other- 
wise flexible backbone. Rigid side chains induced larger 
local fluctuations of the backbone as compared to flexible 
side chains 15[ , yet the persistence length of the entire 
complex was shown to grow super-linear with the length 
of the side chains in contrast to the much weaker de- 
pendency in case of flexible side chains. In other words, 
attaching a large number of rigid side chains yields an 
efhcient way to stiffen a flexible backbone. 

Here, we investigate the conformation of a flexible, 
end-grafted self-avoiding chain of length N with a sin- 
gle rigid side chain (length S) by means of Monte Carlo 
simulations [l^. We find that attaching the side chain 
leads to a considerable anisotropic swelling of the flexible 
backbone and a straightening to a more brush-like con- 
figuration. This phenomenon not only depends on TV and 
S, but also on the position 6 at which the side chain is 
attached. In particular, we show that for fixed N and h 
the squared radius of gyration rises sigmoidally from 
the unperturbed value i?Q with increasing S and levels 
off to a constant for 3> 6/2. In fact, the side-chain 
dependent swelling of the backbone follows a scaling re- 
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monomer b, i.e. by introducing a new type of move. In 
every attempt only those moves are accepted that lead to 
self-avoiding conformations while respecting the impene- 
trability of the substrate. Yet, 'during' a rotation process 
of the side chain the self-avoidance can be violated. In 
connection with the non-ergodicity of all iV-conserving 
algorithms for SAWs [3, this may lead to trapped con- 
formations in the described model. To avoid an ovcrrep- 
resentation of trapped conformations, the simulation is 
continued with a newly generated random conformation 
when the system appears to have run into a dead end. In 
practice, this case occurred for rather short chains only. 



FIG. 1: Two-dimensional sketch of the self-avoiding polymer 
on a lattice as described in the main text. A rigid side chain 
(grey) is attached to the fully flexible end-grafted backbone at 
monomer b ('branching site'). The substrate is impenetrable, 
i.e. the chain conformations are restricted to the half-space 
z > 0. Note that a chain of length N consists of A' + 1 
monomers. 



lation with maximum value A ~ R^b^ /N'^. Contrary to 
an undisturbed end-grafted chain, the ratio ^ond/^g 
creases with S, thus highlighting the successive breaking 
of isotropy, that is also mirrored in dramatic changes of 
the ratios of the principle axes of gyration. Furthermore, 
for b N and S ^ 5/2, the distribution of normalized 
angles (p {{(p) =1) between the surface normal and the 
longest axis of the backbone's gyration ellipsoid follows 
a WeibuU distribution for all combinations A^, b, S. 



II. SIMULATION METHOD 

To simulate a flexible end-grafted polymer chain of 
length A^ (i.e. having A^-|- 1 monomers), we utilize a sim- 
ple cubic lattice with unity lattice spacing (which we have 
taken as our unit of length) . Each monomer occupies a 
single lattice site and self-avoidance is guaranteed by pro- 
hibiting multiple occupations. The polymer is considered 
in a good athermal solvent, i.e. no interactions between 
the monomers besides the self-avoidance are taken into 
account. The monomer with index is placed in the 
plane z = 0, and the halfspace 2 < is chosen as impen- 
etrable while the lattice is taken large enough to avoid 
influences of any other boundary. A rigid side chain of 
length S is attached at monomer b of the flexible chain 
(counted from the grafted end). Figure [1] shows a two- 
dimensional sketch of the model. 

The flexible chain is simulated using a Verdier- 
Stockmayer-type algorithm allowing for kink and 
90°/180°-crank shaft moves (see, e.g., Due to its 

rigidity, the only possible moves of the side chain are ro- 
tations by ±90° and 180° in the xy-, yz-, or a;z-plane 
respectively around the branching site b. A side chain 
rotation can be induced by a regular move (kink, crank) 
involving monomer b or randomly without translocating 



III. RESULTS AND DISCUSSION 

To investigate the behavior of the described polymer 
system, we vary the lengths of the flexible backbone (A^) 
and the rigid side chain (S) as well as the branching site b 
at which the side chain is attached to the backbone. As 
basic read-outs, we then monitor the squared radius of 
gyration Rg, the principle moments of gyration and the 
angle (j) between the longest axis of the gyration ellipsoid 
and the surface normal, defined by cos(j) = V3 ■ e^/L^ 
(cf. Eq. 11])). All these quantities are calculated from 
the backbone monomers only, in order to study the con- 
formational differences of an end-grafted polymer with 
an attached side chain to one without. 

We first determine the squared radius of gyration Rq 
and the orientation angle of the gyration ellipsoid for an 
end-grafted flexible backbone ( A" = 20, . . . , 200) without 
side chain. In agreement with [ll| . we observe an increase 
i?2 ^ ]\f2^ ^i^Yi V K, 0.59 (Fig. Hi). In contrast, the aver- 
age angle of orientation, (i^), only depends very weakly 
on A^ and tends towards an asymptotic value (0) = 55° 
(Fig. [2)d, inset). The entire probability density function 
p((/)) of the orientation angles for A^ = 20 and A^ = 200 
is shown in Fig. [2}d. Clearly, in both cases the mean (0) 
does not coincide with the plateau of the most probable 
value of near to 90°, i.e. the gyration ellipsoid is most 
likely oriented parallel to the substrate ('mushroom'). It 
is noteworthy, that the plateau-like behavior of p((/>) only 
emerges properly for large A^, while a too short polymer 
shows a slight decrease of p((/)) for ^ — > 90°. 

When attaching a rigid side chain of length S at po- 
sition b to the flexible backbone (length A^), we observe 
an increase in the radius of gyration that depends on S*, 
b and N . A representative example (A^ = 100) is shown 
in Fig. [3^. The squared radius of gyration increases sig- 
moidally from the unperturbed value i?o with increas- 
ing S whereas the limiting plateau strongly depends on 
b. The gross shape of the curves can be rationalized by 
considering the limiting cases: for 5* — » and 6 — > the 
unperturbed backbone has to be recovered as the side 
chain vanishes or appears to be 'glued' to the substrate 
surface, respectively. For large S any further increase 
of the side chain is not 'felt' by the mushroom- like back- 
bone, i.e. the curve should level off to a constant. In fact. 
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FIG. 2: (a) The radius of gyration increases like Rq ~ N'^''' 
~ 0.59) for a flexible end-grafted chain (without side 
chain), (b) The probability density function p(0) of the ori- 
entation angle of the polymer's gyration ellipsoid shows a 
plateau-like behavior for (j) —> 90° in the limit of large A''. 
Inset; The average angle {(j)), i.e. the first moment ofp(<l)), de- 
pends only weakly on A*' with an asymptotic value (0) = 55° . 




0.6 5/^0.8 



12 
10 
8 
6 
4 
2 



-4 o 



o 
V o 



■ >< 

r 



+ 



■ . 1_ 



0.0 



0.2 



_l 1 1 1 L_ 



0.4 



0.6 



0.85/^^1.0 



FIG. 3: (a) The squared radius of gyration (Rg) increases 
sigmoidally with the length of the side chain S/N and levels 
off at increasing values with increasing b (here: A'" = 100). 
(b) All curves collapse to a single master curve when apply- 
ing the scaling Eq. (I2J. Open black symbols: b/N = 0.25, 
N = 20,60,80,100,120; filled grey symbols: b/N = 0.5, 
N = 20, 50, 80, 100, 120; crosses: b/N = 0.75, A^ = 20, 100. 



we are able to collapse all data for the side-chain induced 
difference A = i?^ — i?Q to a single master curve (Fig.jSh) 
for various combinations of N and b by the scaling: 

As a result of Eq. ([2]), we clearly observe that A ~ 
Rlb^/N"^ in the plateau region. Thus, for b ^ N (side 
chain attached to the free end) the gain in length fol- 
lows approximately A ~ iV^ which is reminiscent of the 
behavior of flexible end-grafted chains in the 'regime of 
stretched chains' observed in polymer brushes 0]. 

In Fig. m we have plotted the ratio between the mean 
squared end-to-end distance and the mean squared ra- 
dius of gyration, 9 = (^ond)/(^s)' ^ function of S. 
For a simple, free self-avoiding walk this ratio can be cal- 
culated to be = 6 in the limit — > 00 while 9 w 7.6 
for a chain attached to a solid substrate (see offset in 
Fig. m for S = 0). Attaching a side chain results in a 
further increase of 9, indicating that grows faster 

than (Rg)- Again, the effect is strongest for large val- 
ues of S and b. This observation yields further evidence 
that the side chain, due to imposing an excluded-volume 



constraint, causes the felxible backbone to swell and dis- 
entangle to a more straight conformation. This picture is 
confirmed and refined by the subsequent results obtained 
from the study of the gyration tensor. 

While measures the spherical size of a polymer its 
actual shape is better described by the dimensions of the 
gyration ellipsoid Li, L2, L^- The ratio (L^) : (L|) : (L|) 
resembles the asphcricity of the polymer. Our results 
show that for growing S (L§) : (L^) increases up to a 
plateau at 5 « 0.5A^, that depends on b - analogous to 
For N = 100 and b = 0.57V we obtain a maxi- 
mum ratio of approximately 26 : 1, indicating a dramati- 
cally altered conformation in comparison to a simple end- 
grafted chain, where 14.9 : 1 is found [ll|. The ratio is 
even higher for larger values of 6, that is, the backbone 
takes on a pronounced rod-like shape due to the attached 
rigid side chain. 

The ratio : (Lf) between the two shorter axes 

of gyration shows a different behavior: For b ^ G.5A^ it 
changes only marginally with increasing S and takes on 
a value of about 3.1 : 1 independent of b. In contrast, for 
6 = 0.5A^ a small but appreciable increase occurs leading 
to a saturation value of 4 : 1 for A^ = 100. 
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FIG. 4: The ratio Q = {Rcnd)/{Rl) oi the squared radius of 
gyration and the squared end-to-end vector of the backbone 
increases with increasing length of the rigid side chain S/N 
(here: N = 100). 



We next consider the orientation of the backbone in 
terms of the the distribution of orientation angles p{4>) 
between the longest principle axis of gyration, V3, and 
the surface normal. As can be seen in Fig. [5^ for the 
representative example N = 100, elongating the side 
chain leads to a considerable decrease of the average an- 
gle (0) , indicatcing a more brush- like configuration of the 
backbone. This straightening becomes more pronounced 
when the branching site b is moved towards the free end 
of the backbone, i.e. for large S and b N the back- 
bone deviates from the surface normal by less than 20°. 
Concomitant to the decrease of the average orientation 
angle, the entire distribution changes and assumes 
a more compact, i.e. narrow, shape around the mean 
{(j>) for large S (Fig. [SJd). The decrease in width of the 
distribution accompanying the decrease of the average 
highlights the brush-like conformation of the backbone. 

Indeed, the compact shape of p{4>) is a generic feature 
of the backbone for 5 ^ and large S. This is reflected 
by the fact that all distributions can be collapsed to a 
single curve master curve 



„k-l 



f ^ fey' 
pm cxp 



(3) 



with k w 2, A « 1, when considering the normalized 
angle ip ~ 0/(0) and fixing b/N ~ const and S/N = 
const (Fig.[S]). 



IV. CONCLUSIONS 

We conducted Monte Carlo simulations using a flexible 
polymer backbone, end-grafted to a solid substrate, with 
a rigid side chain attached to it. From the behavior of 
the backbone's radius of gyration and the length and 
orientation of its longest principle axis of gyration, wc 
are able to conclude that attaching a rigid side chain 
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FIG. 5: (a) The average angle (0) between the longest axis 
of gyration and the substrate normal decreases for increasing 
lengths of the side chain S/N, indicating a straightening of 
the backbone. Moving the branching site b towards the free 
end of the backbone strongly enhances this straightening, (b) 
The entire distribution p{(f>) shifts and becomes more narrow 
as S/N is decreased (here: N = 100). 



leads to a straightening of the backbone to a more brush- 
like configuration. Depending on the side chain length 
S and the branching site b, the radius of gyration and 
the gyration component perpendicular to the substrate 
is enhanced while the average orientation tends towards 
the surface normal. The effects of the side chain are 
strongest in the case of large b and S. 

For an undisturbed end-grafted chain (Fig. [2]d) the 
probability density p{4>) shows a plateau for angles near 
to 90°, resembling a coiled conformation (mushroom). 
Attaching a stiff side chain leads to a shift to smaller 
angles and a narrowing of p{4>) (Fig. [Hd). It is found 
that p{(j)) can be scaled using a Weibull distribution. Ac- 
cording to the changes of p(0) the arithmetic mean (0) 
changes starting at about 55° in the case of an undis- 
turbed chain to smaller values when a side chain is at- 
tached (Fig. Eli). 

The behavior of the radius of gyration as well as of the 
ratio (-Rg)/(-Rend) indicates a swelling and straightening 
of the backbone due to the presence of the stiff side chain. 
The backbone's change from a mushroom to a more rod- 
like conformation is also reflected in the ratios between 
the axes of gyration: While the two smaller principal 
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heuristic scaling with maximum value /N'^ which impli- 
cates a more brush-like growth of Rg with the backbone 
length N. 

It will be interesting to examine the influence of a semi- 
flexibility of the side chain to obtain a more realistic rep- 
resentation of the above described HA-aggrecan system. 
To approach the biological setting, e.g. the protective 
extracellular matrix, the phase diagram of end-grafted 
polymers with an attached side chain at varying surface 
densities will be of great interest. Work along these lines 
is currently underway. 



FIG. 6: Normalizing the average angle, i.e. setting (p — (f>/ {^(f>), 
for b ~ N leads to a collapse of all distributions p{^p) for 
various N,S with S/N = 0.3 (symbols). The curve is best 
described by Eq. @ (full line) . 

axes only show minor changes, the longest principal axis 
is increased manifold, e.g. 26-fold for b/N — 1/2. Fur- 
thermore, the increase in the radius of gyration follows a 
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